Exact Results for a Kondo Problem in One Dimensional t-J Model 
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We propose an integrable Kondo problem in a one-dimensional (ID) t — J model. With the open 
boundary condition of the wave functions at the impurity sites, the model can be exactly solved 
via Bethe ansatz for a class of Jr,l (Kondo coupling constants) and Vl,r (impurity potentials) 
parametrized by a single parameter c. The integrable value of Jl,r runs from negative infinity to 
positive infinity, which allows us to study both the ferromagnetic Kondo problem and antiferromag- 
netic Kondo problem in a strongly correlated electron system. Generally, there is a residual entropy 
for the ground state, which indicates a typical non-Fermi liquid behavior. 
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Recently, considerable attention has been drawn by the theory of impurities in both the Fermi and Luttinger 
^ \ liquids, and many new developments have been addressed. This renewed interest in the quantum impurity problems 
was partially stimulated by the search for non-Fermi-liquid fixed points beyond the well known Luttinger liquid and 
CZ) its plausible relevance to the theory of heavy fermions and cuprate superconductors. A great progress is the study on 
the multi-channel Kondo problem, which has been extensively studied by the conformal field theory in the presence of 
boundary |l],|2j and Bethe ansatz ^-|) , These powerful methods allow to obtain the low temperature thermodynamics 
near the critical point. Using simpler bosonization and renormalization group techniques, Kane and Fisher have 
| shown that a potential scatter center embedded in a Luttinger liquid is driven to a strong-coupling fixed point by 
£^ ■ the repulsive electron-electron interactions 01. This is the first time to show that a single impurity in a Luttinger 
liquid behaves rather different from that in a Fermi liquid, and directly stimulates the study on the problem of local 
perturbations to a Luttinger liquid and especially on the Kondo problem in a Luttinger liquid. The Luttinger-Kondo 
problem was first considered by Lee and Toner JjJ , who found the crossover of the Kondo temperature from power 
law dependence on the Kondo coupling constant to an exponential one. Subsequently, a poor man's scaling was 
performed by Furusaki and Nagaosa who addressed a conjecture which states that ferromagnetic Kondo screening 
may occur in a Luttinger liquid. The boundary conformal field theory || gave out a classification of critical behavior 
for a Luttinger liquid coupled to a magnetic impurity (without impurity potential). It turns out that there are only 
^sO • two possibilities, a local Fermi liquid with standard low-temperature thermodynamics or a non- Fermi liquid observed 
by Furusaki and Nagaosa. It is now clear that the non-Fermi-liquid behavior is induced by the tunneling effect of 
conduction electrons through the impurity which depends only on the bulk properties but not on the detail of the 
impurity 

Despite these important progress, the problem of few impurities (potential, magnetic, especially both) embedded 
in a strongly correlated system is still not well understood. We remark that there are a few progress related to this 
I ■ problem: A spin S > 1/2 impurity in a spin 1/2 Heisenberg chain solved many years ago by Andrei and Johannesson 
and generalized to arbitrary spins by Lee and Schlottmann, and Schlottmann jl3| ; an integrable impurity in the 
supersymmetric t — J model |14| with a very complicated hamiltonian solved by Bedurftig et al. We remark here that 
these models are somewhat special for the absence of backward scattering off the impurity, and therefore are hard 
• • ■ to give a deep understanding for the Kondo problem in a Luttinger liquid, where the backward scattering is crucial 
. £h ! t° the fixed point of the system. In addition, Wang and Voit has proposed an integrable model of a single magnetic 
. impurity in a 5— potential Fermi gas [jllj with a special value of Kondo coupling constant. 

Attempting to understand effects of magnetic impurities in a strongly correlated electron system, we study the 
properties of the integrable t — J 

model coupled to two magnetic impurities sited at the ends of the system. Our starting point is the following 
hamiltonian 

H = H o + Hi , 

N a -l N a -l 

H = -J2 ( C l C J+i° + h - c -) + 2 E Pi ' S ->+i + Vn jnj+ i], (I) 

J=1,CT j = l 

Hi = JlSi • S L + Vlui + JR,S Na ■ S R + Vnn Na , 
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where Cj a {C^ a ) are annihilation (creation) operators of the conduction electrons; V, Jr,l, Vr,z are the near- 
est neighbor interaction constant, the Kondo coupling constants and the impurity potentials respectively; Sj = 
\ Eiru 1 Cj a <j aa 'Cj a ' is the spin operator of the conduction electrons; Sl,r are the local moments with spin-1/2 sited 
at the left and right end of the system respectively; rij — C^Cji + C^Cji are the number operator of conduction 
electrons; N a is the length or site number of the system. Notice that the single occupation condition rij < 1 is 
assumed for the hamiltonian (1). We remark that the model is very resonable for the absence of redundant terms in 
the hamiltonian. 

It is well known that Ho is exactly solvable for V = —1/4, 3/4 |l5| , ^6| . In this letter, we only study V — 3/4 case 
while both the charge- and spin-sector can be described by a Luttinger liquid. The V — — 1/4 (supersymmctric t — J 
model) can be followed without any difficulty. By including the impurities, any electron impinging on the boundaries 
will be completely reflected and suffer a reflection matrix Rj l or Rj,R. 01 • The waves are therefore reflected at either 
end as 

e ik jX ^ R j L (kj)e- ik t x , x ~ 1, 
e «y* _» R~ l R (k j )e- ik i*-* ik i N *, x~N a (2) 

The reflecting Yang-Baxter equation |l7] , |l8f| 

Sji(q 3 - Qi)Rj,a(Qj)Sji(qj + qi)Ri, a (qi) = Rt,a(qi)Sji(qj + qi)Rj,a{<lj) s jMj ~ Qi)i a = R,L, (3) 

constrains the integrablity of the present model. Here Sji is the electron-electron scattering matrix in the bulk. By 
evaluating eq.(3), we obtain that the present model is exactly solvable with the following parametrized J a and V a 
(a = L,R): J a = — 8/ (2c a + 3)(2c a — 1), V a — (4c^ — 7)/(2c a + 3)(2c a — 1), where Cl,r are two arbitrary real constants. 
Generally, Jl and Jr may take different values. Here we consider only Jr = Jl, Vr = Vl, i.e., cl = cr = c case. The 
solution of our model in the integrable line is quite similar to those of other integrable models with open boundaries 
0,0. Notice that the reflection matrix Rl.r (K± in ref. [17]) is an operator one rather than a c- number matrix. 
The spectra of the hamiltonian (1) are uniquely determined by the following Bethe ansatz equations 

Qj + i Qd- *( c - !) r i±i 9j -rqi+i ^ qj - rX a - | 

X a -i(c- |) A a +i(c+ |) 2 -j-j- yr Kj2 rq] + I = TT TT a « - r\ p + i 

X X a + i{c-\)X a -i{c+\y iyi^-rft-i iii/i^-^-*' 

with the eigenvalue of (1) given by E — 2N — J2j=i 4/(4q| + 1). Here iV is the number of conduction electrons, 
qj = l/2tan(fc J /2), and kj and A a are the rapidities of charge and spin respectively. M is the number of spin-down 
electrons. 

Below we discuss the ground state properties for different regions of parameter c. 

(i) c > 1. The system falls into the ferromagnetic Kondo coupling regime and no bound state can exist at low 
energy scales. The ground state is thus described by two sets of real parameters {qj} and {A Q }. Define the quantities 

M N 

z c NM) = -{- Q M + ^f\<t>M- E *i(?-Aa)+ E fcfo -©)]}> 

a a=-M j=-N 

N M 

^ a (A) = — {0,(A)~ E i( A -*)+ E ^(A-A Q )}, (5) 

with the phase shifts c (g) = -6 2 (q) +4tan- 1 [ (? /(c- 1)], ^(A) = -6» 2 (A) + 4tan- 1 [A/(c+ 1/2)] -4tan- 1 [A/(c- 1/2)] 
induced by the impurity in the charge- and spin-sector respectively, and n (x) = — 2 tan _1 (2a;/n). Note above we have 
used the reflection symmetry of the Bethe ansatz equations to include solutions with q_j = —qj and A_ Q = — X a . The 
Bethe ansatz equations are solved by Z^^qj) — Ij/N a and Zff (A) = J a /N a , where Ij and J a are non-zero integers. 
In the ground state, Ij and J a must be consecutive numbers around zero symmetrically to minimize the energy. The 
roots qj and A a becomes dense in the thermodynamic limit, and we define their densities as 

c , s MnM s m dZ kW f({ . 
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The cutoffs of q and A in the ground state are ±Q and ±A respectively, which correspond to Z C N (±Q) = ±(N + 
l/2)/N a , Z s Na (±A) = ±{M + 1/2) /N a . In the thermodynamic limit N a -> oo, N -> oo, and N/N a finite, we find 
that the energy is minimized at A — > oo, which gives a result of N = 2M by integrating eq.(6). This leaves a spin 
triplet ground state which is contradicting to the Furusaki-Nagaosa conjecture ||. Notice that the magnetization is 
given by 1/2 (AT+2-2M). 

(ii) 1/2 < c < 1. The system falls also into the ferromagnetic Kondo coupling regime. However, unlike case (i), two 
bound states of electrons can be formed around the impurities in the ground state, which correspond to two imaginary 
q modes at q = ±i(l — c). The Bethe ansatz equations for the real modes are thus reduced to 



~ i ^Ng = r gj + i( y c - l ) 1j - ic Qj + *( c ~ 2 ) 
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qj + | qj - i(c - 1) q 3 + ic qj - i(c - 2) 



x 



r=±l Ijtj q l Tqi + 1 o=l ^ - rA *-2 

\ a -i(c- f) \ a +i{c+\) 2 yr X a - r qj + § = „ „ A a - r\ fj + z 
Aq, + i(c — |) Aq. — i(c + |) r i£ /i A Q - r % - § j£ 11 A Q - rA^ - i ' 

Following the same procedure discussed above, we obtain again a spin triplet ground state. This can be understood 
in the following picture: The bounded electrons and the local moments form two spin-1 local composites due to 
the effective attraction and the ferromagnetic Kondo coupling. However, the itinerant electrons impinging on these 
composites will screen their moments partially due to the indirect Kondo coupling induced by the electron-electron 
correlation. In such a sense, we recover Furusaki-Nagaosa's conjecture |||, though the local moments are not completely 
screened. 

(iii) —1/2 < c < 1/2. The system falls into the antiferromagnetic Kondo coupling regime. No bound state appears 
in the ground state. By integrating eq.(6) in the thermodynamic limit, we have N + 2 = 2M. This indicates a spin 
singlet ground state, a similar result to that of the Kondo problem in a Fermi liquid. 

(iv) — 1 < c < —1/2. No bound state exists in the ground state. It seems that the ground state should be a spin 
triplet. We note that both J a and V a take positive values in this region. The repulsive boundary potential dominates 
over the Kondo coupling. It will repel the conduction electrons to form singlet with the local moments and thus make 
the Kondo coupling ineffective. In such a sense, no Kondo screening occurs, which strongly indicates that both the 
boundary potential and the electron-electron correlation in the bulk have significant effects to the Kondo problem in 
a Luttinger liquid. 

(v) —3/2 < c < —1. For this case, the system is still in the regime of antiferromagnetic Kondo coupling but with 
a weaker or attractive boundary potential. Two imaginary A mode at A = ±i(c+ 1/2) and two imaginary q modes 
at q = ±i(c + 1) appear in the ground state. These modes correspond to the formation of bound singlet pairs of two 
conduction electrons with the local moments. By integrating the densities of the real modes in the thermodynamic 
limit, we still arrive at a spin singlet ground state. 

(vi) c < —3/2. The Kondo coupling is ferromagnetic and the boundary potential is repulsive. There is no bound 
state in the ground state. Following the same procedure discussed in (i), we obtain again N = 2M for the ground 
state. Therefore, there is no Kondo screening in this region, which contradicts to the Furusaki-Nagaosa conjecture. 

The thermodynamics of the present model can be calculated in a closed form based on the Bethe ansatz equations 
(4). This allows us to obtain the temperature and magnetic field dependence of the free energy which contains three 
parts of contributions, i.e., the bulk term, the boundary term and the Kondo effect term. Here we omit the detail 
of calculation following the standard procedure which can be found in some excellent works P,|l9|-pT|. The Kondo 
effect induced free energy is the most interesting one which takes the following form at low temperatures (Hereafter 
we assume c is an integer or half integer) 

^ Hl + C| Bl |(A)]dA /• ln[l + C|n 2 |(A)]dA 

- -T^n( ni ) J 2cQsh(7rA) ~ Ts W n{n 2 ) J 2cosh(?rA) , (») 

where n\ and 72.2 are two c-dependent integers, Cn(A) are elements of the following coupled integral equations 

In 77(A) = £o(A ]T - - ([1]G + [2]) ln[l + Gln[l + Ci(A)], 

lnC„(A) = G{ln[l + Cn+i(A)] + In[l + C«-i(A)], n > 1, 
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lnCi(A) = -Gln[l + rr 1 (A)]+Gln[l + C 2 (A)], (9) 
lim {[n] ln[l + C„+i(A)] - [n + 1] ln[l + C„(A)]} = — , 

n — >oo J 

where fi is the chemical potential and £o(A) = 2 — 4A 2 /(4A 2 + 1); [n] and G are integral operators with the kernels 

(A) = (n/2-7r)/[A 2 + (n/2) 2 ] and l/[2 cosh(7rA)] respectively; H is the external magnetic field. Notice that we have 
omitted the excitations breaking the bound states in deriving (9) for the energy gaps associated with them. Their 
contributions to the free energy are exponentially small at low enough temperatures. 

For H = and T — > 0, the driving term Gln[l + r] _1 (A)] in (9) is nothing but T _1 e s (A) with e s (A) the dressed energy 
of the spin waves which has the asymptotic form e s (A) = 2-7reo exp(— 7r|A|) for |A| — ► oo, where eo is the energy density 
of the ground state. This allows us to formulate the low-temperature expansion of (9). The asymptotic solution of 
(9) is given by functions Cn(x), x = ln[(7reo) /T] — tt\X\ which are monotonically decreasing in x for all n and tending 
to finite limits £ rl+ as x — > oo. The limits are given by Cn+ = sinh 2 (nH/T)/ sinh 2 (H/T) — 1. No anomaly appears in 
the free energy of the bulk and the open boundary. They contain only constant terms and T 2 terms up to the order 
0(T 2 ). An interesting feature is that with different c value, the system may show both Fermi- and non- Fermi-liquid 
behavior. To show this, we study the entropy of the ground state for a variety of c regions. For c > 1 or c < —3/2, 
there is no bound state in the ground state. n\ and ni take the values 2|c| + 1 and 1 — 2\c\ respectively. This leaves 
a residual entropy of the ground state as S g = ln[(2|c| + l)/(2|c| — 1)]. For c = 0, the local moments are completely 
screened. Both n\ and ri2 are equal to unit. Thus the entropy of the ground state is zero and the system flows to 
a Fermi-liquid fixed point. For c = —1/2, n\ = 2 and ri2 — 0. The residual entropy takes a value of ln2. While for 
c = —1, rii = —1/2, ri2 = 3/2, S g = ln3. Notice that for c = —1/2 and c = — 1, the system falls into the regime of 
antiferromagnetic Kondo coupling but with a non-zero residual entropy. We remark that c = 1/2 and c = —3/2 are 
two critical points since at these points, both J a and V a are divergent. The residual entropy has different limits for 
c — » 1/2 + + (— 3/2 + + ) and c — > 1/2 + 0~(— 3/2 + CP). From the above discussion we conclude that the system 
generally has a c-dependent residual entropy which strongly indicates a non-Fermi-liquid behavior. The system can 
flows to a Fermi-liquid fixed point only in a narrow parameter region of c ~ 0. Such a fascinating effect can be 
understood in the following picture: The charge-spin coupling induced by the backward scattering off the impurity 
introduces an effective boundary field to the local moment, which means the charge degrees of freedom join the Kondo 
effect. However, unlike a real magnetic field, the "effective field" does not split the degeneracy of different orientation 
of the local moment. In fact, when the conduction electrons impinging and leaving the impurity, the incident waves 
and the reflection waves feel different strengths of the impurity spin. One is 1/2 — c and the other is 1/2 + c or vise 
versa. The finite residual entropy is a result of cooperative effect of the charge- and spin sectors. 

In summary, we introduce an integrable model of Kondo problem in a ID strongly correlated electron system. With 
different values of the parameter c, the system can show either a Fermi liquid behavior or a non- Fermi liquid behavior 
beyond that obtained by Furusaki and Nagaosa Q . A non-singlet ground state in a antiferromagnetic Kondo coupling 
region is obtained. It is found that the residual entropy depends not only on the self magnetization of the ground 
state, but also on the interaction parameter c, which we interpret as a "cooperative effect" of the Kondo coupling and 
the impurity scattering. It is instructive to apply renormalization group analysis and conformal field theory to reveal 
a fullcn picture for such an interesting problem. 

Acknowledgment 

One of the authors (YW.) is supported in part by the National Natural Science Foundation of China and the 
Natural Science Foundation for Young Scientists, Chinese academy of Sciences. 



[1] J.L. Cardy, Nucl. Phys. B 240, 514(1984); B 324, 581(1984);D. Lewellen and J.L. Cardy, Phys. Lett. B 259, 274(1991). 

[2] A.W.W. Ludwig and I. Affleck, Nucl. Phys. B 428, 545(1994); I. Affleck in Springer Series in Solid State Sciences, 118, 

ed. by A. Okiji and N. Kawakami (Springer- Verlag 1994) and references therein. 

[3] N. Andrei and C. Destri, Phys. Rev. Lett. 52, 364(1984). 

[4] A.M. Tsvelick and P.B. Wiegmann, J. Stat. Phys. 38, 125(1985). 

[5] N. Andrei and A. Jerez, Phys. Rev. Lett. 74, 4507(1995). 

[6] C.L. Kane and M.P.A. Fisher, Phys. Rev. Lett. 68, 1220(1992); Phys. Rev. B 46, 15233(1992). 

[7] D.-H. Lee and J. Toner, Phys. Rev. Lett. 69, 3378(1992). 

[8] A. Furusaki and N. Nagaosa, Phys. Rev. Lett. 72, 892(1994). 



4 



[9] P. Frojdh and H. Johannesson, Phys. Rev. Lett. 75, 300(1995). 

[10] Y. Wang, J. Voit and F.-C. Pu, Phys. Rev. B 54, 8491(1996). 

[11] Y. Wang and J. Voit, Phys. Rev. Lett. 77, 4934(1996). 

[12] N. Andrei and H. Johannesson, Phys. Lett. 100 A, 108(1984). 

[13] K.-J.-B. Lee and P. Schlottmann, Phys. Rev. B 37, 379(1987); P. Schlottmann, J. Phys.: Cond. Matt. 3, 6619(1991). 

[14] G. Bediirftig, F.H.L. E/31er and H. Frahm, Phys. Rev. Lett. 77, 5098(1996). 

[15] P. Schlottmann, Phys. Rev. B 36, 5177(1987). 

[16] P.A. Bares, G. Blatter and M. Ogata, Phys. Rev. B 44, 130(1991). 

[17] E.K. Sklyanin, J. Phys. A: Gen. Math. 21, 2375(1988); C. Destri and H.J. de Vega, Nucl. Phys. B 361, 361(1992); B 374, 
692(1992). 

[18] A. Foerster and M. Karowski, Nucl. Phys. B 396, 611(1993); B 408, 512(1993). 

[19] C.N. Yang and CP. Yang, J. Math. Phys. 10, 1115(1969). 

[20] M. Takahashi, Prog. Theor. Phys. 46, 1388(1971). 

[21] N. Andrei, K. Furuya and J. H. Lowenstein, Rev. Mod. Phys. 55, 331(1983). 



5 



